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C^l ' Abstract 

Gamma-spectrometric methods suitable for the characterization of highly enriched 
'~. , uranium samples encountered in illicit trafficking of nuclear materials are presented. In 

■ particular, procedures for determining the ^^*U, ^^^U, ^^*U, ^^^U and ^^'^U contents 



and the age of highly enriched uranium are described. Consequently, the total uranium 
. content and isotopic composition can be calculated. For determining the ^^®U and ^^^U 

contents a low background chamber was used. In addition, age dating of uranium was 
also performed using low-background spectrometry. 
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^ ■ 1 Introduction 



The detection of illicit trafficking of nuclear materials, including "nuclear smuggling" is nec- 
essarily followed by forensic investigations in order to determine the possible origin of the 
Q ■ material. These investigations involve, beside classical forensic analysis, also the techniques 

^ \ of nuclear forensic science. Nuclear forensic investigations are carried out by laboratories 

C ' capable of handling and analyzing nuclear material of unknown origin, or "nuclear forensic 

laboratories" . 

At least three levels of nuclear forensic capabilities can be distinguished ^ |2] ■ These 
are: categorization, characterization, and full attribution of materials to the point of loss-of- 



■ control 



Categorization means the determination of the type of the investigated material ("low- 
enriched" uranium, "high-enriched" uranium, plutonium, irradiated materials, other radioac- 
tive substances), in order that a threat assessment can be done. Categorization is usually 
performed by high-resolution gamma-spectroscopy, possibly using portable equipment so that 
a fast threat assessment can be done on the spot where illicit trafficking is detected. 

Characterization involves a more precise and more detailed determination of isotopic 
composition, crystal structure, trace-element analysis, age dating, particle size analysis etc. 
Some of the methods used for characterization of nuclear material of unknown origin are 
mass spectroscopy, alpha spectroscopy, X-ray diffraction analysis, electron microscopy, etc. 
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Upon characterizing nuclear material of unknown origin, some questions related to the 
production processes and intended end use may be addressed. If an appropriate knowledge 
base is also available, full attribution of materials to the point of loss-of-control might be 
possible. 

In this work we wish to show that simple gamma spectroscopy can be used not merely 
for categorization, but also for basic characterization of highly enriched uranium-'^ (HEU) of 
unknown origin. In particular, it is shown how gamma spectroscopy can be used for very 
accurate determination of the age and the isotopic composition of HEU. These are properties 
which are usually assessed by mass spectroscopy, which is relatively complicated, requires 
sophisticated (and expensive) equipment and special sample preparation is needed. On the 
contrary, gamma spectroscopy is relatively simple and, since it is non-destructive, it does 
not require any special sample preparation. Therefore, although gamma-spectroscopy is not 
capable of all types of analysis done by mass spectroscopy (e.g. trace element analysis), 
in nuclear forensics it might be the preferred choice for laboratories that do not possess an 
appropriate mass spectrometer or do not need the additional information supplied by mass 
spectrometry. 

In 2001 the International Technical Working Group for Combating Illicit Trafficking of 
Nuclear Materials (ITWG) organized an inter-laboratory comparison exercise (Round Robin 
Exercise - RRE) in order to assess and compare means and methods of forensic analysis on 
samples of seized nuclear material In the Exercise samples containing about 2 g of 

a highly enriched uranium oxide powder were sent to the participating laboratories. The 
participants were asked to determine the properties of this material, relevant for its forensic 
analysis. The results were expected to be returned to the organizers in periods of 24 hours, 
1 week and 2 months. 

In this paper the gamma-spectrometric methods which can be used to determine the 
234u, 235u^ 238u^ 232u^ and ^^^\] , as wcll as the ^^^Th and 24iAm contents and the age of 
HEU are reviewed. In order to determine the age of the sample, the activity of ^^^Bi was 
measured. The applicability of the methods is demonstrated by determining these quantities 
for the sample received within the Round Robin Exercise. 

Three high-purity germanium (HPGe) detectors were used for the measurements. The 
^^''U and 235^ content of the sample were determined using a large-area planar high-purity 
germanium detector. A procedure for self-absorption correction was applied and the count 
rates per unit mass of the relevant gamma lines of ^■^^U and ^ssy -^gj-g determined. The 
isotope contents were then calculated by comparing these count rates to the corresponding 
ones obtained from the measurement of a reference sample. 

A coaxial HPGe detector in a low background iron chamber was used to measure the 
activities of ^^^U, ^'^^U, ^•^^Th and ^^''Bi, based on the knowledge of the absolute efficiency 
curve of the detector in the relevant energy interval. 

Finally, the activities of ^sey ^nd ^''^Am were estimated from the spectra taken by a 
medium-area planar HPGe detector, using an intrinsic efficiency calibration method. 

From the above measurements the activity ratio ^^'^Bi/^'^^U was calculated, and using this 
result the age of the sample was also determined. The presence of ^^ay ^nd 236y indicates 
that the investigated material was produced by reprocessing spent nuclear fuel. 

At the end of the paper a discussion of the results is given, together with a comparison 
of the measured values to those obtained using mass-spectrometry by other laboratories 
participating in the Round-Robin exercise. 

^In this work the term "enriched uranium" is used to denote uranium with a relative isotope abundance 
in 235 u larger than the natural 235 y content of uranium. The expression "highly enriched uranium" (HEU) 
has various meanings in the literature, but here we will only use it to denote isotopically enriched uranium 
with more than 80% 235u, contrary to most of its uses elsewhere. 
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2 Determination of the ^^^U and 235^ content 

The 234u and ^a^U contents of the Round Robin Material ("RRM") were determined using; 
a 20 cm^ planar HPGe detector (Canberra GL2020 with active diameter of 50.5 mm and 
thickness of 20 mm). A reference material ("RFM") was used for calibration. The RFM 
was U3O8 powder of known isotopic composition (^s^U: 90.6 ± 1.5%; ^^^V: 8.35 ± 0.20% 
and 23'*U: 1.02 ± 0.07%). This means that in 1 g of the RFM there is a total of 0.846 g of 
uranium consisting of 0.766 ±0.12 g ^ssu, 0.070 ±0.001 g ^^^V and 0.0086 ± 0.0004 g of ^a^u. 
The reference material was transported to Hungary before 1960, so its age at the time of the 
measurement (June 2001) was more than 40 years. Both nuclear materials (the RRM and 
the RFM) were placed separately into a cylindrical polyethylene container of 2.9 cm inner 
diameter. 

The spectra of the samples were measured at a 10.8 cm distance from the detector. Four 
measurements with amounts of approximately 0.5, 1, 1.5 and 2 g were carried out with each 
material. The count rates of the gamma-transitions at 120.9 keV (^^''U) and 185.7 keV (^^'^U) 
were standardized in unit mass as the ratio, K(m,E), of the count rate to the total used 
mass of the sample, to, and then they were plotted versus m (see Fig. 
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Figure 1: Count rate of 120.9 and 185.7 keV per unit mass of the sample as a function of 
the total mass of the sample. 



The count rates per unit mass at energy E, K(m, E), were fitted with the fimction 



Kim,E)^Ko{E)- 



aETTi 



(1) 
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which models the law of self-absorption within the sample. The parameter Ko{E) is the 
value of the curve at ?7i = and physically it represents the net count rate at energy E of 
1 g sample, corrected for self-absorption. The parameter as provides information about the 
matrix of the sample. The ^'^^U and ^ssy contents of 1 g of the sample were derived by 
comparing the count rate per unit mass of the RRM, Ko{E), to that of the RFM, Kq{E). 
More precisely, the ^^"^U content of the Round-Robin material, M{U23A), was calculated by 
inserting the corresponding data into the formula 

«<--'^-'<--'§^' ,2, 

where the primed quantities correspond to the reference material. An analogous formula was 
used for the the ^35^ content. A non-linear least squares fit was performed to obtain the 
parameters Koil20.9keV), Klj{120.9keV), KQ{185.7keV), and i^^(182.7fcey). In this way 
the contents of ^a^u and ^asy in 1 g of the RRM were found to be 0.00807±0005 g/g and 
0.754±0.011 g/g, respectively. 



3 Investigating the sample by low-background gamma- 
spectroscopy 

3.1 Determination of the ^38^ content 

The samples of the RFM (1.05 g) and the RRM (0.975 g) were held in a closed cylindrical 
polyethylene container of 2.9 cm inner diameter as in the measurements described above 
and assayed at a distance of 6.5 cm from a 150 cm'^ coaxial high-purity germanium detector 
(PGT PIGC 3520 with 34.1 % relative efficiency) in a low background chamber. Their 
gamma-spectrum is shown in Fig. Altogether with the background spectrum. 

The gamma-peaks at 766.37 and 1001 keV (Fig. ^1 and Fig. ^2) of 234"ipa can be 
observed clearly. Since ^'^^'"Pa is a daughter of ^^^U, its 1001 keV gamma-transition was used 
for estimating the ^ss^j content. The count rate per unit mass at 1001 keV corresponding 
to the RRM was compared to the one corresponding to the RFM. Because the sample is 
"thin" and practically transparent for the 1001 keV line, there was no need to correct for 
self-absorption. Therefore, the count rates per unit mass at 1001 keV Ko{1001keV) and 
KQ{1001keV) were evaluated by simply dividing the measured count rate by the total mass 
of the sample. Inserting the received values (see Table^ into a formula analogous to formula 
(O, the ^^^U content of 1 g sample was estimated to be 0.065±0.002 g/g. 



3.2 Age determination 

As reported in Refs. |H1E1) the activity ratio ^^^Bi/^'^^U can be used as a chronometer for 
uranium age dating. The gamma-transitions at 609.3, 1120.3 and 1764.5 keV were identified 
for ^^"'Bi (Fig. |21;l-c3). The activity of ^^"'Bi in Ig of the sample is calculated by the formula: 

where K is the count rate of Ig of the sample, B is the emission probability and e is the 
detector efficiency measured by point-like standard sources. Because of using a thin sample, 
the count rate does not need to be corrected for self-absorption in the sample. 

The activity of ^^''Bi in the RRM was calculated in Ref. ||5| and it is summarized in Table 
13 ^^^Bi is a daughter of ^'^"'U, which decays through ^'^"Th to ^^^Ra, which, in turn, decays 
to ^^''Bi through three short-lived nuclides. If the activity of ^^''Bi, A, and the activity of 
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Figure 2: Gamma-spectrum of the background (1), the reference material (2) and the Round- 
Robin material (3) 
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234 



U, Aq, are known, the age T of the sample can be derived from the foUowmg formula 



1 



^Th^RaT 



(4) 



A 

A 2' 

where Xrh = 0.288 x lO^^^g-i ^^d A^q = 13.86 x lO^^^g-i ^re the decay constants of ^^^Th 
and 226Ra |3, respectively. The age of the RRM was estimated in this way to be 23 ± 3 
years. The age of the RFM was also calculated and it was found to be 42 ± 3 years, which 
was the expected value. 



3.3 Determining the activity of ^^^U and '^^'^Th 

During the operation of a nuclear reactor the isotope 232y formed. Therefore, if the 
presence of 232]j jg detected in the investigated material, it is probably a re-processed material 
produced from spent nuclear fuel. The activity of ^ssy jg determined from the activity of its 
gamma-emitting daughters, ^^^Bi and ^"^Tl. However, since these isotopes are also present 
in the decay chain of ^"^^Th (Fig. |2J), determining the activity of ^32^ ^^so involves the 
determination of the activity of ^^^Th. 
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Figure 3: Decay scheme of ^'^^Th and ^'^^U. 



Observing the decay scheme of ^szy and 232^^]^ (Fig. 13), and the gamma energies emitted 
by their daughters, it can be seen that the activity of ^^^Th can be calculated from the 
activity of ^^^Ac, which, in turn, can be determined from the gamma peaks of ^^^Ac at 
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911.204 and 968.971 keV. Using the law of radioactive decay one obtains 

Aac228 = ARa228 = ^Th232(l " e-^«"^^«*) (5) 

where Aac22?,i ^_R,a228 and ATh232 denote the corresponding activities, XRa228 is the decay 
constant of ^^^Ra, while t is the age of the sample. 

The activity of ^32^ ^g^j^ calculated from the activity of ^^^Bi and ^°^T1, which can be 
determined from the count rates at their gamma peaks. The presence of these isotopes in the 
sample, however, may due both to the decay of ^^^Th and to the decay of ^^^\J. Therefore, 
the activity of ^^^Bi and of ^os^j ^g^^^ ]-,e calculated as the sum of two terms, one of them 
accounting for the buildup from ^'^^Th and the other term accounting for the buildup from 

232^. 



ABi2l2 — Z — ^Th232 



A 



(7232 



^ _ ATfe228e'^^°^^"* - XR.a228e-^^ 
A/?a228 ^ At/i228 



1 — At/232/ ATh228 



(6) 



where / = 35.9% is the decay branching probability of the ^^^Bi 208rpj decay (see Fig. 
Using the measured activities of ^^^Ac, ^^^Bi and ^''^Tl, the measured age of the sample and 
the known half-lifes, the activities of ^ssy g^j^j 232 rpj^ ^-^^ obtained by solving equations lO 
and © for Au232 and ATh2Z2- 

The activities of ^-'^^Bi and ^osrpj -^gre determined from the measured count rates using 
the known efficiency of the detector (in the same way as the activity of ^^'^Bi in section 
and they are presented in Table |2] The activity of ^^^Ac can be estimated from its peaks 
at 911.204 and 968.971 keV. After subtraction of the background, however, the net count 
rates of the 911.204 and 968.971 keV peaks turned out to be zero within the measurement 
error, both in the spectrum of the RRM and of the RFM. Therefore, the upper limits of 
these count rates may be taken to be equal to their statistical errors. This leads to an upper 
limit for the value of the '^^'^Th. activity of about 0.37 Bq/g in the RRM and 0.19 Bq/g in 
the RFM. Inserting these values into the first term in formula ((HJ, one obtains 0.40 Bq/g and 
0.19 Bq/g as the contribution of the ^'^^Th decay to the activity of ^^^Bi in the RRM and 
in the RFM, respectively. Since these values are much smaller than the uncertainty of the 
measured ^^^Bi and ^os^p^ activities, the first term in equation © can be safely neglected. 
Therefore, inserting the weighted average of the measured values of ABi2i2 and Arnm/ f 
into the left-hand side of formula lO one obtains the ^32^ activity as 419 ± 8 Bq/g for the 
RRM and 42 ± 1 Bq/g for the RFM. 

4 Determination of ^36^ content and ^^^Am activity by 
a medium area planar HPGe detector 

The spectrum of the RRM was also taken by a medium area planar HPGe detector, ORTEC 
SGD-36550-S, with a crystal diameter of 36.7 mm. The 49.37 keV peak observed in the 
spectrum was identified for ^^^U decay (Fig. Although this peak is very close to the 
49.550 keV peak of ^^*U, in the present case this does not disturb the evaluation of the 
49.37 keV peak. Namely, in the observed sample the isotopic ratio of ^38^ jg gQ thai 
the intensity of its 49.550 keV peak is less than the background level, i.e. it is below the 
detection limit. Therefore, the 238u peak at 49.550 keV can be neglected when evaluating 
the one at 49.37 keV. 
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The ratio of 236u g^^^j 235jj contents was evaluated by comparing the yield of the 236jj 
peak at 49.37 keV to the yield curve obtained from the yields of the 58.6, 84.2 and 90 keV 
peaks of ^^^U and it was found to be (0.67 ± 0.10) x lO'^. Then the ^ssy content was 
calculated to be (0.51 ± 0.08) x 10"^ g/g for the RRM. 




-| 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

40 45 50 55 60 65 70 

Energy (keV) 

Figure 4: Gamma-spectrum of the Round-Robin material taken by a medium-area planar 
HPGe detector. 

Similarly, the 59.6 keV line was identified for ^''^Am and its activity was calculated by 
peak ratio technique as above and found to be 125 ± 15 Bq/g for the RRM. 

In the spectra of the RFM the peaks of ^36^ g^j^^j ^^^Am could not be identified. More 
precisely, the upper limit of the count rates at the peaks of ^ssy g^Yid ^'''^ Am may be estimated 
as the statistical error of the count rate in the areas of the spectrum where these peaks should 
be (at 49.37 keV and at 59.6 keV). In this way, for the RFM the upper Hmit of the ^^'^U 
content was estimated to be about 0.0008 g/g, while the ^^^Am activity in the RFM was 
estimated to be not more than 15 Bq/g. 

5 Discussion 

In the present paper independent NDA methods for HEU characterization are described. 
These methods were developed within the two-months time frame of the Round Robin Ex- 
ercise organized in 2001 by the International Technical Working Group for Combating Illicit 
Trafficking of Nuclear Materials OS]. In fact, the results for the 235jj content and total 
uranium content were included already in the 24-hours report. The results of the detailed 
calculations of the ^^*U, ^'^^U, ^'^^U, 238^ and 232^ contents were used for preparing the final, 
two-months report, together with the results of mass spectrometric measurements. The age 
of the Round-Robin material determined by the method described here was also available 
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before the end of the two-months time frame of the exercise. The resuhs of the presented 
research are compared to the results obtained by mass spectrometry in Table |3| 

Age dating is one of the most interesting results of the presented research. In the Round 
Robin Exercise in 2001 the age of HEU was estimated by gamma spectrometry for the first 
time. This method does not require the use of calibration standards and its uncertainty is 
comparable to that of mass-spectrometry (see Table 121. Age dating requires measuring the 
content of ^'^''U and ■^^'^Bi as precisely as possible. Unlike commercially available software 
OEIE] (see also ^^), the procedure described in this work can provide a quantitative 
estimation of the ^^*U content which is sufficiently precise for age dating. 

The estimation of total uranium content is an important result of this analysis. The 
obtained value of 0.835 ± 0.015 g/g is close to the uranium content of UaOs (0.846 g/g). 
But if we note that the amounts of the three uranium isotopes (^^"'U, ^ssy g^y^^ 238y-j ^-j^^ 
RRM are significantly smaller than in the REM (Tabled or Eig. Qand Fig. Elbl-b2), which 
is known to be UaOg, it follows that the matrix of the RRM is different from the matrix of 
the REM and it cannot be identified as UaOg . 

The ability to determine a small ^ssy amount in a low background chamber makes it 
possible to precisely determine the isotopic composition of a material even with nearly 100% 

of 235U. 

Estimating the 232jj gj^^ 236u content provides important information about unknown 
nuclear material. Because the nuclides of ^32^ ^nd ^se^j gj,^ produced during the operation 
of a nuclear reactor, their existence in the investigated sample shows that the HEU material 
was produced from reprocessed spent nuclear fuel. A good energy resolution planar HP- 
Germanium (small or medium area) should be used to observe the 49.37 keV peak of ^36^ 
for improving the accuracy of the result. 

The accuracy of determining the age of HEU and the ratio of major uranium isotopes 
j-234u^ 235JJ g^j^^ 238^-^ gamma spectroscopy is comparable to the accuracy of mass spec- 
troscopy. Eurthermore, gamma-spectroscopy can also provide information on the total ura- 
nium content of the sample, which cannot be determined by mas spectrometry. The isotopic 
ratio of ^'^^U, however, which is easily measured by mass spectrometry, is difficult to resolve 
by gamma spectroscopy, and below a certain limit (around 0.1%) its gamma-spectroscopic 
determination becomes impossible. On the other hand, only gamma spectroscopy can be 
used to measure small activities of ^^^U. Since ^36^ and 232]j are produced during the oper- 
ation of a nuclear reactor, identifying the presence of either of them is sufficient to conclude 
that the sample was produced from reprocessed spent nuclear fuel. 

Ideally, for the full attribution of nuclear materials of unknown origin a balanced inter- 
action between non-destructive and destructive techniques should be achieved. It has been 
shown here that high-resolution gamma-spectroscopy is a powerful tool of non-destructive 
analysis for the characterization of highly enriched uranium. Since gamma-spectroscopy is 
one of the most simple and most cost effective techniques, it might be the first choice for 
the basic characterization of highly enriched uranium. If a more detailed characterization, 
including, e.g., full trace element analysis, description of crystal structure, estimation of par- 
ticle size distribution etc. is needed, other methods of assay, if available, should be called 
for. 
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Tables 



Table 1: Uranium content and isotopic composition of the Reference Material and of the Round- Robin Material. (* Activity (Bq/g), * uranium 
content of UsOg.) 



Nucleus 


Reference Material 






Round-Robin Material 




Isotope ratio (%) 


Content (g/g) 


K (cps/g) 


Ko (cps/g) 


Content (g/g) 


Isotope ratio (%) 




1.02(4) 


0.0086(3) 


4.9(2) 


4.6(2) 


0.0081(4) 


0.98(7) 


235 u 


89.5(1.0) 


0.757(12) 


163.0(7) 


160.5(7) 


0.745(15) 


89.8(7) 


238 u 


9.5(2) 


0.080(1) 


0.00351(3) 


0.00312(3) 


0.0715(35) 


8.6(5) 


236 u 


< 0.1 


< 0.0008 






0.0051(8) 


0.6(1) 


232](j 




5.4(2) X 10-11 


42(1)# 


419(8)# 


5.4(1) X IQ-i" 




Total 




0.846* 






0.830(20) 





Table 2: Activity of ^■^^U, ^^"^Bi, ^^^Bi and ^"^Tl. Emission probabilities were taken from Ref. [7|, except for 120.9 keV of ^'^'^U, which was taken 
fromRef. 0. 



Round-Robin Material Reference Material 



Nucleus 


Energy 
(keV) 


Emission 
probability 

(%) 


Detector 
efficiency 

(%) 


Count rate 
per unit mass 
(xlOO cps/g) 


Activity per 
unit mass 
(Bq/g) 


Count rate 
per unit mass 
(xlOO cps/g) 


Activity per 
unit mass 
(Bq/g) 




120.9 


0.0342 


0.730(15) 


460(15) 


1.84(6) X 10^ 5 


490(15) 


1.91(6) X 10" ig 




609.3, 1120.3, 1764.5 






Averagc= 


2.00(25) 5 


Avcragc= 


7.06(15) 5 




727.3 


6.58 


0.483(1) 


13.6(1) 


429(15) 


1.38(5) 


43(3) 




583.0 
860.3 
2614.3 


84.5 
12.14 
99 


0.58(1) 
0.45(1) 
0.15(2) 


76.2(2) 
8.37(8) 
25.5(1) 

Average= 


156(4) 
153(6) 
172(25) 
155(4) 


7.72(15) 
0.76(5) 
2.46(5) 
Average= 


16(1) 
14(1) 
17(3) 
15(1) 



Table 3: Isotopic composition, 232y content and the age (in 2001) of the Round-Robin Material, measured by different laboratories using 
mass-spectrometry (MS) and gamma-spectrometry (GS) [30]. 



Laboratory nickname 




235^ 


256^ 


238 "u 


■^•i^U (g/g) Age (years) 


Method 


Azores 


0.97 


89.99 


0.68 


8.37 


22.2-22.6 


MS,GS 
MS 


Barbados 




85.6 ±3.8 








GS 


Borneo 


0.85 ± 0.15 


86.7 ± 1.5 


0.57 ± 0.08 


11.9 ± 0.9 




MS 


Chatham 


0.960 ± 0.001 


89.94 ± 0.06 


0.0643 ± 0.003 


8.462 ± 0.006 




MS 


Galapagos 


0.96 


89.89 


0.68 


8.47 




MS 


Mindanao 


0.96 ± 0.40 


89.9 ± 0.11 


0.678 ± 0.23 


8.443 ± 1.29 


22.4 ± 1.2 


MS 


Tobago 


1.05±0.07 


89.37±1.8 


0.69±0.05 


8.88±0.2 




MS 


Tonga 


0.967± 0.001 


89.99 ± 0.02 


0.679 ±0.001 


8.362 ± 0.005 




MS 


Trinidad 


0.955±0.075 


90.01±0.35 


0.673±0.030 


8.365±0.033 


23.5± 0.5 


MS 



Present work 0.98±0.07 89.8±0.7 0.6±0.1 8.6±0.5 (5.4 ±0.1) x lO^^" 23±3 GS 



